Connection preserving deformations and g-semi-classical orthogonal polynomials 



Christopher M. Ormerod, N. S. Witte and Peter J. Forrester 

Abstract. We present a framework for the study of q-difference equations satisfied by q-semi-classical orthogonal 
systems. As an example, we identify the g-difference equation satisfied by a deformed version of the little q-Jacobi 
polynomials as a gauge transformation of a special case of the associated linear problem for g-Pyi- We obtain a 
parameterization of the associated linear problem in terms of orthogonal polynomial variables and find the relation 
between this parameterization and that of Jimbo and Sakai. 
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pH | 1. Introduction 

Monodromy representations have been a central element in the study of integrable systems [5]. A pioneering 
step was the parameterization of the condition that a linear second order differential equation with four regular 

CN ■ 

■ singularities {0,t, l,oo} has monodromy independent of t, in terms of the sixth Painleve equation Pyi 

o 



This theory was elaborated upon by Garnier 125| and Schlesinger |59j . and culminated in the 1980's with the studies 



and 18\. For a matrix linear differential equation of the form 



of the Kyoto School 1371 1381 139) . A contemporary perspective of the theory can be found in the monographs 

o 

(1.1) ^-Y(x)=A(x)Y(x), 
' where 



A(s)=£— ' 

L — ' X — Ot{ 

i 

one expects the general solution to be multivalued with branch points located at a = {oi}. By evaluating a solution 
on any element of the homotopy classes of closed loops, [7], in some manner around a selection of the poles, one 
obtains the equation 

F( 7 (l))=y( 7 (0))M [T] . 

This relates solutions on different sheets of a Riemann surface. The set {Mpyi : 7 : [0, 1] — > C} is a representation of 

the fundamental group of the compliment of the poles, r = 7Ti (CPi \ {a}). The aim is to deform the linear system 

with respect to a chosen deformation parameter, t, so that the representation of Y does not depend on t. In the 

theory of monodromy preserving deformations, a natural choice of parameters are the poles of A. This leads to the 
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classical Schlesinger's equations [59 
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Given a 2 x 2 linear system with four poles, {0, i, l,oo}, and zero of Ai%(x) at y, imposing the isomonodromic 
property in the variable t requires that y satisfies Pyi |21j . 

The relationship of the theory of monodromy preserving deformations to orthogonal polynomials arises because 
under certain conditions on the orthogonality measure the polynomials and their associated functions form an 
isomonodromic system, albeit one with a particular restriction. Under fairly general conditions the derivative of 
each polynomial in the system is expressible in terms of linear combinations of other members of the orthogonal 
polynomial system, an observation first made by Laguerre [45] . The conditions for when this is the case have been 
given in some generality by Bonan and Clark [12] , and by Bauldry [7J. In particular, a semi-classical weight will 
satisfy these conditions and the three term recurrence tells us that we may express the derivative of the polynomial 
in the system of orthogonal polynomials as a rational linear combination of the polynomial itself and the previous 
polynomial in the system. The rationality of this linear problem means that such orthogonal polynomial systems 
satisfy a linear problem of the form Ql.ip . 

The notion of a semi-classical weight or linear functional was introduced by Maroni [50] as an attempt to 
characterize the classical orthogonal polynomials in a coherent framework and guide the quest of looking for systems 
beyond this class. By appropriately extending the work of Laguerre, Magnus |49j was able to show that a semi- 
classical, deformed orthogonal polynomial system parameterized a special case of the monodromy preserving 
deformation considered by Fuchs [21] . This allows one to express special solutions of Pyi in terms of coefficients 
of orthogonal polynomial systems. Conversely, this also allows key quantities relating to orthogonal polynomials to 
be expressed in terms of solutions of Pyi- In addition to the various determinantal solutions of integrable systems 
provided by the theory of orthogonal polynomials, the application of integrable systems to orthogonal polynomials 
have resulted in advances in the calculation of various statistics of interest in random matrix theory (see e.g. |19j ). 

For g-difference equations, an analogue of the theory of monodromy preserving deformations is the theory of 
connection preserving deformations 40j . The linear problem of interest is given by the n x n matrix equation 



Instead of the basic information being contained in the relation between the value of solutions on different sheets of 
the Riemann surface, the variables of interest are associated with the relation between two fundamental solutions, Yq 
and , which are holomorphic functions at and oo respectively. Much of the theory concerning the existence of 
these solutions has remained relatively unchanged since the pioneering days of Birkhoff and his followers [U I10L 114] . 
If these solutions exist, then one may meromorphically continue these solutions on C, and furthermore form the 



(1.2) 



Y(qx) = A{x)Y(x), 



where 



A(x) = A + A x x + ... + A m x 
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connection matrix, P{x), specified by 



Y (x) = F 00 (x)P(x), 



which is quasi-periodic in x |10j . From the Galois theory of q- difference equations, which primarily considers a 
classification of problems of the form (|1.2|) . we know that the entries of P(x) are expressible in terms of elliptic 
theta functions |58J. 

In the same manner as monodromy preserving deformations, one may consider a deformation of (|1 .2[) that 
preserves the connection matrix. An appropriate choice of deformation parameter turns out to be the roots of the 
determinant of A and the eigenvalues of Aq and A m . By considering a 2 x 2 linear system with m — 2 and choosing 
the deformation parameter to be proportional to two of the roots of the determinant and the two eigenvalues of 
Aq, Jimbo and Sakai [40 showed that the connection preserving deformation was equivalent to a second order 
g-difference equation admitting the sixth Painleve equation as a continuum limit. 

We have remarked that semi-classical orthogonal polynomial systems give rise to monodromy preserving de- 
formations relating to Painleve equations. A natural problem then is to investigate the relationship between q- 
scmi-classical orthogonal polynomial systems, connection preserving deformations, and the q-Painleve equations. A 
number of different approaches to constructing isomonodromic analogues for the difference, g-difference and elliptic 
equations of the Sakai Scheme |57j have been proposed recently pH |4j 1551 165] . which differ in varying degrees 
from what we offer here. One other work which is close to the spirit of the present work is that of Biane [9] . 
However there is a history of studies into g-semi-classical orthogonal polynomial systems which was not motivated 
by the above considerations. Shortly after the introduction of the semi-classical concepts Magnus extended this 
to the g-difference systems and in fact to the most general type of divided difference operators on non-uniform 
lattices in a pioneering study |48j. In addition Maroni and his co-workers have made extensions to difference and 
g-difference systems in a series of works |43L 1511 1421 1521 127j . These later authors have successfully reproduced 
parts of the classical Askey Tableaux (which was achieved most fully by Magnus at the level of the Askey- Wilson 
polynomials) however the application of their theoretical tools beyond the classical cases have invariably been 
made to specialised or degenerate cases and failed to make contact with the discrete and g-Painleve equations. 
A slightly different methodology has been the approach of Ismail and collaborators |33L 1311 1161 132) . who have 
derived difference and g-difference equations for orthogonal polynomials with respect to weights more general than 
the semi-classical class, much in the spirit of the Bonan and Clark and Bauldry studies, and so the matrix A{x) 
is no longer rational. This approach has not been applied to systems beyond the classical Askey Tableaux, and 
consequently not made contact with the discrete Painleve systems. The most recent work of Biane |9|, and of 
Van Assche and co-workers [641 111] has addressed some of the shortcomings discussed above, however while these 
authors have uncovered the spectral structures of the theory they have yet to elucidate the deformation structures 
required. It is our intention to complete this task by laying out the deformation structures. 

Our contributions in this paper are to first formulate an extension of the classical work of Laguerre for finding 
differential equations satisfied by orthogonal polynomials, when the differential operator is the g-difference 



(1.3) 



Dq, Xi f(xi,..., 
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This is done in Section 4, after first having introduced preliminary material from orthogonal polynomial theory, and 
connection preserving deformations in Sections 2 and 3 respectively. In Section 5 we apply this extension to the 
specific case of a deformation of the little q- Jacobi polynomials [2] . We give a parameterization of the associated 
linear problem in terms of variables relating to the orthogonal polynomial system. However, this contains redundant 
variables, and in fact a set of three natural coordinates can be identified which suffice to parameterize the linear 
problem. When written in terms of the natural coordinates, the linear problem implies the q-Pyi equations 

, n - a 7 a 8 (y2 ~ ait)(y 2 - a 2 t) 

(l-4a) 2/1 yi =- 



(1.4b) 2/22/2 



(2/2 -a 3 )(y 2 -o 4 ) 
Q3«4(2}i - «5*)(2/i - ast) 
(2/1 ~a r ) (2/1 - as) 



where 2/j = Vi(t) and j/j = yi(qt). We show that this has the consequence of implying the r-functions have determi- 
nantal solutions in terms of Hankel determinants of the moments of the little q- Jacobi weight. Also we show that 
the three term recurrence, written in terms of the natural coordinates, manifests itself as a Backhand transformation 
which relate to a translational component of the extended affine Weyl group of type D^. Throughout we shall 
assume that q is a fixed complex number such that < |g| < 1. 



2. Orthogonal polynomials 

Our starting point is a sequence of moments, {pk}kLo- F rom this we define a linear functional, L, on the space 
of polynomials, where L(x ) = Hk- An orthogonal polynomial system is a sequence of polynomials, {p n }^Lo, such 
that p m is a polynomial of exact degree m and 

(2.1) L(piPj) = Sij. 

In other words, these polynomials are orthonormal with respect to the given linear functional. This condition 
defines the coefficients of p n for all n so long as the Hankel determinants consisting of the moments fio, ■ ■ ■ , A*2n, 
given in (HQ)) , do not vanish [TH [61] . 

In the case of the classical continuous orthogonal polynomials, this linear functional, L, is typically some integral 
of the multiplication of the argument with some weight function over some support. Linear functionals associated 
with discrete orthogonal polynomials are specified by a weighted sum, such as Jackson's g-integral |28L I36|, 162), 163] . 
Any orthogonal polynomial system where (|2.1[) holds satisfies the classical three term recurrence relation, given by 

(2.2) a„ + ip n+ i = {x — b n )p n - a n p n -i- 
We parameterize the coefficients of these polynomials by 

Pn{x) = ~f n X n + "fn.lX n ~ 1 + "fn,2X n ~ 2 + • • • + Jn.n- 
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It is possible to determine all the coefficients, and hence the a n and b n , in terms of the jitfc's [171 161] . We set 



(2.3) 



A„ = det 



for n > 1, with Aq = I, and 



(2.4) 



S„ = det 



/ Mo Mi 

Ml M2 

Mn-2 Mn-1 

\Mn-l Mn 



Mo Mi 
Mi M2 



Mn-l \ 
M» 



M2n-2/ 



M«-2 Mn 
Mn— 1 M"+l 



for n > I , where So = 0. Then we have 

(2.5a) 

(2.5b) 

(2.5c) 



\Mn-l Mn ■•■ M2n-3 M2n-1/ 



A n _|_i 

A„_i A n+ i 
= A 1 ' 

A n +i A ra 



as given in [171 161j . 

Given a sequence of valid moments, one may define the Stieltjes function 



/ = E 



f^nX 



-rt-l 



ri=0 



We define the associated polynomials and associated functions by the formula 



fp 



-1 + £nj 



where 4> n -i is a polynomial and e„ is the remainder. The orthogonality condition implies that e„ ~ 7 n a; ™ 1 + 
0(x~ n ~ 2 ). In fact, by using (|2.2I) , it is possible to find the large a; expansions for these polynomials in terms of the 
a n and b ni giving 



n-l 



' n — 2 n — 1 n — 1 



(2.6a) 



(2.6b) 



=-„ i a?- - /- 1 X! 6 < + x "~ 2 E E ^ - E ° 2 ] + °(^ 3 ) ] . 
1 x-*- 1 + 2 5>+ 3 ( E E ^ + E a * ] + oc*-- 



i=0 



i=0 j'=0 



i=l 



where in the second equation use has been made of the large x expansion of j n e n = (x — b n ) 1 7„-ie n -i + [x — 
bn)^ 1 0>n+i1n+i e n+i [30L I17L 1611 149] . Utilising this we can write some explicit relations between the coefficients, 
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7„ and 7 n ,fc's, and the a n and b n , 

(2.7a) a,, - 7 "" 1 



(2.7b) b 



•n-l 



In 

7«-l,l 7n,l 



7n-l 7n 

which hold for n > 1. By equating (fp n )Pn-i with (fp n -i)Pn we have the relation 

(2.8) 0n-lPn-l - 4>n~2Pn = ^n-lPn - ^nPn-X = , 

a n 

which is polynomial by the left hand side and where the final equality follows from (|2.6[) . 

It is clear that the sequence of functions, {e n }^ =0 , i s a solution of (|2.2|) that is independent of {p n }- We find 
it convenient to introduce the matrix 



(2.9) Y n 



Pn e n /w 

The three term recursion relation is equivalent to the relation 



(2.10) Y n+1 = M n Y n 

where 

M n = 



({x- b n ) 



\ 



n+l O-n+1 

1 



3. Connection preserving deformations 

In this section we revise the established classical theory of systems of linear g-difference equations [lj 114] . 
The general theory concerns the m x m matrix system 

(3.1) Y(qx) = A(x)Y(x), 

where A(x) is rational in x. We call A the coefficient matrix of the linear q-difference equation. One may easily 
verify that such an equation possesses two symbolic solutions, namely, the infinite products 

(3.2a) A(x/q)A(x/q 2 )A(x/q 3 ) 
(3.2b) A(x)- 1 A(xq)- 1 A{xq 2 )- 1 . . . , 

which do not converge in general. We may suitably transform the problem so that A(x) is polynomial, which we 
parameterize by writing 

A(x) = A a + A x x + . . . + A n x n . 

The matrices Ao and A n are assumed to be semisimple with eigenvalues pi, ... , p m and Ai, . . . , A m respectively. 
Regarding the solutions of (|3.1[) . we present the following theorem due to Carmichael 

Theorem 3.1. Suppose the eigenvalues of Aq and A n satisfy the condition 
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then there exists two solutions, called the fundamental solutions, of the form 

Y =%x D °, 

Yqq YqqQ 2 X , 

where Yq and are holomorphic functions in neighbourhoods of x — and x = oo respectively and Dq and Dqo 
are diag(log 9 pi) and diag(log g Aj) respectively and u — log ? x. 

We may use (|3.1[) to continue both solutions meromorphically over C \ 0. From these solutions, we define the 
connection matrix to be 

(3.3) P(x) = Y OQ (x)- 1 Y (x). 
The evolution of P9x) in x is given by 

P(qx) =Y oo {qx)- 1 Y (qx), 

=Y oo (x)- 1 A(x)- 1 A(x)Y (x), 

(3.4) =P{x)- 

Hence this function is q-periodic in x. 

These fundamental solutions may be related to (|3.2p via a conjugation of transformations of (|3.1[) such that the 
solutions given by (13.21) converge. Hence (|3.2[) gives us information regarding the roots of the determinant of the 
connection matrix. If the Zi are the zeros of det A(x), then Y^ and Yq are possibly singular at {q n+1 Zi : n G N} 
and {q~ n Zi : n € N} respectively. Therefore we expect the poles and zeros of the determinant of P(x) to be g-power 
multiples of the Z{. 

In the situation of monodromy preserving deformations we introduce a parameter, t, into A and consider 
what conditions on the evolution of t are required so that the monodromy representation is preserved. It was the 
innovation of Jimbo and Sakai |40j to introduce a parameter t in a manner that preserves the connection matrix, 
P (x), through the evolution t — > qt. By the observation noted in the above paragraph regarding the poles and 
roots of the determinant of P(x), we may infer that the latter are preserved in the shift Zj —¥ qzj. However, in 
doing this, we must also consider the eigenvalues of Aq and A n to be parameters. This suggests that the zeros of 
the determinant of A and eigenvalues of Aq and A n are appropriate choices of parameter for the connection matrix 
preserving deformation. 

We parameterize pa = det A by letting a subset of the roots be constant in t, while the other roots are simply 
proportional to t. By supposing P(x, t) = P{x, qt) we arrive at the implication 

P(x,qt) =Y oo (x,qt)- 1 Y (x,qt) ^ Y^x^t^Yoix^t) = P{x,t), 

which defines a matrix, _B, via 



Yoix^tjYoix^)- 1 =Y QO (x,qt)Y 00 (x,t)- 1 = B(x,t). 
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Since both Yo and Yqq are independent solutions this leads to the necessary condition that the evolution of Y must 
be governed by a second linear equation, 

(3.5) Y(x,qt) = B(x,t)Y{x,t). 

Conversely, it is easy to see that if Y satisfies (|3.5j) . then the evolution t — > qt defines a connection preserving 
deformation. 

Since Y satisfies an equation in x and an equation in t this imposes the necessary compatibility condition 

(3.6) A(x, qt)B{x, t) = B{qx, t)A{x, t). 

This compatibility condition implies a (/-difference equation satisfied by pb — deti?, 

/ .s PA{x,qt) 

PB(qx,t) = — - — —p B {x,t), 

PA{X,t) 

which may be solved up to a factor of a function of t. We shall assume that ps is rational in x. Other informa- 
tion regarding asymptotic behavior of A gives us specific information regarding the form of B. Furthermore the 
compatibility condition and the determinantal constraints often results in an overdetermined system allowing us 
to construct a representation of B in terms of the entries of A and hence, find (/-difference equations in t for the 
entries of A. 

We do not pursue this line explicitly. Rather, in the following sections we will show how the (/-difference 
equations satisfied by a particular deformed g-semiclassical orthogonal polynomial system leads linear systems 
satisfying (|3.1[) and (|3.5[) . 

4. (/-difference equations satisfied by orthogonal polynomials 

4.1. The ^-difference calculus and (/-special functions. A reference for the (/-difference calculus and also 
the /i-differential calculus is a book by Kac |41) . We first recall some of the basic properties in relation to (/-difference 
equations. The first property is the (/-analog of the product and quotient rule, given by 

(4.1a) D q , x {fg) =jD qtX g + gD q . x f, 

=.fD q , x g + gD qyX f, 

(4.1b) /7\ u'K.f fiK.,<j _ 

\gj 99 



where / = f(x) = f(qx). Associated with the (/-difference operator is the antidcrivative, known as Jackson's 
(/-integral, as defined by Thomae and Jackson |36l 1621 163] . We express the definite integral of Thomae [621 163j by 

(4-2) / f(x)d q x=(l-q)y^q k f(q k ). 

Jo fc=0 

This was subsequently generalized by Jackson |36j to 

pb r-b pa 

(4.3) / f(x)d q x= f(x)d q x- f(x)d q x, 

Ja JO JO 

where 



/ f(x)d q x = a(l-q)Y / q k f(aq k ) 
Jo k=o 
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If / is continuous, then 



Urn / f{t)A q t = / f(t)dt. 
i^Jo Jo 



We introduce the g-Pochhammer symbol 



(a;q) n = (1 - o)(l - aq) . . . (1 - aq^ 1 ), 



and 



We also adopt the notation 



(a;g)oo = (1 - a)(l - aq)... 



(ai,a 2 , . . . , a m ; q) n = (ar, q) n (a 2 ; q) n ■ ■ ■ (a m ;q) n - 



The symbol (a; q)^ is often called the q-exponential as 



-a(aqx;q) c 



(1-9) ' 

We are now able to express Heine's basic hypergeometric function |29j . as it was re- written by Thomae 



(4.4) 



291 



a, b 
c 



q;t 



y- tm (a, b l g)r. 



Also relevant is the integral formula for the basic hypergeometric function 



(4.5) 



a. b 



2<Pl 



q;t = 



(M;«) £ 



1 palL-i (xta,xq;q)oo 

X la &l — — ^ d q X. 



(1 - q)(c,q;q)oo J Q (xt,f;q) 



Jacobi's elliptic multiplicative theta function, as defined by Jacobi's triple product formula, may be expressed 



a.s 



This satisfies the equation 



Of importance is the g-character, 



Qq{z) = ( q,-qz,—;Q 



9 q (qz) = qzO q (z). 



e q (x/ c ) ' 



satisfying 



(4.6) 
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4.2. q-difference equation in x. We shall henceforth assume that the log g-derivative of the weight specifying 
the linear form of the g-orthogonal polynomial system is rational, parameterizing its g-derivative via the equation 

(4.7) W{x)D q>x w{x) = 2V(x)w(x), 

where W(x) and V(x) are polynomials in x. This is the g-analogue of the notion of semi-classical weight functions 
familiar in the theory of classical orthogonal polynomials [50J and which was proposed by Magnus |48j and subse- 
quently by others [51|, 143 j . Henceforth we regard systems satisfying these conditions as q-semi-classical orthogonal 
polynomial systems. 

Lemma 4.1. Assuming w satisfies (|4.7|) . the Stieltjes function, f , satisfies the q-difference equation 

(4.8) W(x)D qtX f(x) = 2V(x)f(x) + U(x), 
where U is polynomial such that degU < degV . 

We call the polynomials W, V and U the spectral data polynomials. We will use (j4.8[) as the basis for our 
derivation of the x-evolution for the g-semi-classical orthogonal systems. Using different methods, this has been 
derived in 1331 H6|, 1314 1321 for general (i.e. beyond semi-classical) g-orthogonal polynomials. 

Theorem 4.2. [481 142) The matrix Y n satisfies 

(4-9) D qx Y n A n Y nj 



where 

(4.10) A 



Q n - V -a„6n 



(W(x) - 2x(l - q)V(x)) [a^.x tt n _i -V-(x- b n ^)Q n ^) ' 
and Q n and Q n are polynomials specified by 

(4.11a) 9„ =W(e n D q , x p n - p n D q>x e n ) + 2Ve n p n , 

(4.11b) Cl n =a n W(e n - 1 D q , x p n - p n -iD q>x e n ) 

Pn—l^n 

) - 2Vx{\ - q)a n e n ^ 1 D q 

,xPn ■ 

PROOF. Using P~gj) and 

„ _ <fin-l £n 
Pn Pn 

we find 

/ K-i 



IT-n I "-'i 2V( Pn-l „ T „„ ( £n\ , 2Ve n 
WD q . x U = -WD q . x — M . 

V Pn J Pn \Pn J Pn 

This may be expanded to give 

WipnDq^n-x - 4> n -iD q , x p n ) - 2U0„_ip„ - Up n p n _ W{e n D q , x p n - p„D q<x e n ) + 2Ve n p ri 



PnPn PnPn 

Setting this equal to Qn/{PnP n ) defines O n as 

(4.12) 0„ =W(p n D qtX <f> n _i - <j) n -iD q . x p n ) - 2V<f) n - 1 p n - Up n p n , 

=W(e n D qtX p n -p n D qtX e n ) + 2Ve n p n . 
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The first expression is a linear combination of polynomials which verifies that 0„ is a polynomial for all n. The 
second expression is one in which every occurrence of a p n or its derivative is multiplied by a factor of e„ or its 
derivative. Hence, by (|2.6j) . there exists an upper bound for the degree of the polynomial 0„ in x that is independent 
of n. 

Using the fact that 0„ is polynomial and 

(4.13) a n 4> n -lPn-l ~ a n <Pn-2Pn = 1) 

we have 

{a>n<t>n-lPn-l ~ a„0„_ 2 Pn)S n = W(j) n D q>x (f) n - 1 ~ 4> n -\^> q,xPn) ~ 2U0„_ip„ - Up n p n . 

By appropriately equating factors divisible by 4> n -\ on one side and factors divisible by p n on the other side, we let 
f2„ be the common factor by writing 

p„(/) n -ifl n =<p n -l (a n 0nPn-l +(W - 2Vx(l - q))D q ^ x p n + Vp n ) , 

=Pn {a n (p n - 2 @n + WD q<x (j) n -i - V(f) n -i - Up n ) . 

The first equality is a rearrangement of the required q-difference equation for p n . The second expression for Q n is 
equivalent to 

(4 14) ft = a «0»-20» + WDg^n-! _ y _ Up n 

4>n-l <j>n-l <f>n-l 

We use the expression for 0„ in terms of n _i and p n to give 

_ Wa n (j) n -2PnDq,x4 l n~l _ W a n 4> n - 2 4>n-lD q/xPn _ 2Va n 4>n-24'n-lP n 



Ua n (/) n -2P„p n t WD q , x <t> n -i Up n 



<fin-l <fin-l 9n-l 

which, by using the equality a„p„0„_2 = a-nPn-i<t>n-i — 1 to eliminate occurrences of 4>n-2, is equivalent to 

fin = a n W(p n ~iD q ^4> n ^i - 4> n ~2D q ,xPn) - U(2a„0„_ 2 p„ + 1) - Ua n p n -ip n . 

This expresses il„ as a linear combination of polynomials and hence Sl„ is a polynomial. To obtain the large x 
expansion, by dividing the first expression for f2„ by 0„_ip„ we have 

e„ + WD q 

,xPn y (1 - q)VD q 

,xPn 



Pn Pn Pn 

Using (|4.11a[) , we find that £l n may be written as 

,xP?i ^n 

Wp n -ip n D q 

,x^n . 

2a n p n -iVe n p n WD q 

,xPn -r r 

2x{\ - q)VD q 

,xPn 

\ L n — ~r i ~r V • 

Pn Pn Pn Pn Pn 

Using f|4. 13|) to cancel the factors of p n -i gives (|4.11bj) . 
A rearrangement of (|4.15l) is 

(W - 2x(l - q)V)D q , x p n = (fi„ - V)p n - a„9„p„_i, 

where fl„ and Q n are given by (|4.11ap and (|4.11bj) . Mapping n — > n — 1 in this relation and expressing p n -2 in 
terms of p n and p n -i gives 

(W - 2x(l - g)U)£) g]X p„_i = a„e„_ip„ + (0„_x - V - (x - fo„-i)6„-i)p„-i, 
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showing that the orthogonal polynomials form a column vector solution of (|4.9[) . 

To see that e n /w satisfies the same (/-difference equation, we consider W D qtX (j> n _i + WD qiX e n — WD QtX (fp n ), 
which we reformulate as 

WD a + WD a =Wp n D q>x f + WfD q 

=2Vfp n + Up n + f ((fi„ - V)p n - a„6„p„_i + 2V Pn - 2Vp n ) , 
= [Up n + (0„ + V)(j) n -x - a„e„0 n _ 2 ] + [(^n + V) e n - a n <d n £n-i} ■ 
Hence by subtracting WD q>x (f) n -i, defined by (|4.14p . we find 

WD qtX e n = (f2 n + V)e n - a„e„e„_i, 

and from this 

wD qiX e n - e n D qiX w 



D q 



WW 

(n„ + v) 2Ve n 
= w w 

w{l-2x{l-q)^f) 

CO V) hi — n f) 6n ~ 1 
= W-2x(l-q)V ' 

while the shift n — > n — 1 and (12.21) gives the compatible evolution equation for e n -\/w. □ 

Since p n is of degree n and the leading order term of e„ about x — oo is x~ n ~ , we immediately obtain upper 
bounds for the degrees of 8 n and Sl r , 



(4.16a) deg9„ < max(deg W - 1, degV - 2,0), 

(4.16b) deg Q n < max(deg W, deg V - 1, 0). 



It follows that knowledge of W and may be used in conjunction with (|2.6p to determine f2 n and 0„ in terms of 
sums and products of the cii and hi. 

We remark that (I4.9[) can be rewritten in a manner more familiar in the context of connection matrices where 
Y n is a solution of the linear g-diffcrcnce equation 

(4.17) Y n {qx) = {I- x(l - q)A n )Y n (x) = A n Y n (x). 

As revised in Section 3, under appropriate conditions on A n established by Carmichael |14j . or the more general 
conditions of Adams pQ, this equation permits two fundamental solutions: Yoo,n and Y"o,n- We form the connection 
matrix by using the fundamental solutions of (I4.9[) in (|3.3I) 

(4.18) P n (x) = (Y^x^Yo^ix). 

Solutions of (|4.9p also satisfy (|2 . 10[) and so the transformation n — > n + 1 is a connection preserving deformation. 
This also implies that the connection matrix is independent of the consecutive polynomials chosen in the orthogonal 
system. Therefore the connection matrix is an invariant of the orthogonal polynomial system. 
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In further specifying the evolution in x for Y n , we have a compatibility relation between the evolution in x and 
the evolution in n. This compatibility relation, defined by equating the two ways of evaluating D q)X Y n +i, namely 
D qtX M n Y n = A n+ iM n Y n , results in the condition 

(4.19) M n {qx)A n (x) + D q , x M n (x) = A n+1 (x)M n (x). 

The entries of the first row are equivalent to the recurrence relations of Magnus [481 142] 

(4.20a) (x - b n )(n n+1 - n n ) =W - x{\ - q)(Q n + v)- a 2 n e n ^ + a 2 n+1 Q n+1 , 

(4.20b) Q n -i - fl n+1 ={x - fo n _i)6„_i - (qx - b n )Q n . 



Using the relation (|2.8[) we find 

(4.21) detY» = ~ p "- l£ " = _L, 

w a n w 

From this we may deduce 

detY n (qx) _ w{x) _ W(x) 



(4.22) dct(7 - x(l - q)A n ) = 



det Y n {x) w{qx) W(x) - 2x{l - q)V(x) ' 
or equivalently 

(4.23) ,(1 - q) det A n - TrA n = w _^_ q) y , 
and these in turn imply the additional recurrence relation 

(4.24) x(l - q) (n 2 n -V 2 - a 2 n e n Q n ^) = {{x - 6„_i)6 n _ 1 - - n n )(w - x{l - q)(Q n + V)). 

The compatibility between the evolution in x and n for more general orthogonal polynomials has given rise to 
associated linear problems for discrete Painleve equations |53L 164] . Many of these associated linear problems are 
differential-difference systems |34j . That is to say that the evolution in x is defined by a differential equation, while 
the evolution of n is discrete. The first occurrence of a discrete Painleve equation in the literature is thought to 
have been deduced in this manner 60j. 

The 2x2 linear problem derived for orthogonal polynomials is one in which the coefficient matrix, A n (x), is 
rational. If we follow the theory of connection matrices, we apply a transformation that relates the linear problem 
in which A n is rational to another linear problem in which the coefficient matrix is polynomial. With respect to 
Birkhoff theory and (|3.1I) . the coefficient matrix obeys the proportionality constraint 

det A oc W(W - 2x(l - q)V). 



That is to say that when referring to the connection matrix for orthogonal polynomial systems, we do not distinguish 
between the roots and the poles of the determinant of the linear problem. 
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4.3. q-difference equation in t. We now turn to the new direction that we wish to present. It is at this 
point we let w(x) = w(x,t), and hence consider polynomials which are both functions of x and t. For functions 
f(x,t) with two independent parameters we will adopt the notation 



f = f(x,t) = f(qx,t), 
f = f&t) = f(x,qt). 

We distinguish a base case in which deg0„ = and deg£!„ = 1, corresponding to degV = 1 and degVF = 2, as 
being completely solvable and a case in which the connection matrix is known [46j . However by suitably adjoining 
^-exponential factors that depend simply on t to the numerator or denominator of w{x) we introduce roots or poles 
into w{qx)/w{x). This has the effect of increasing up the degree of W(x) and V(x). Furthermore it imposes a 
rational character on the logarithmic g-derivative of w with respect to t : 

(4.25) R(x, t)D q>t w(x, t) = 2S(x, t)w(x, t), 

where R(x, t) and S(x, t) are polynomials in x. These cannot be arbitrary polynomials in x as there is an implied 
compatibility condition. This arises because there are two ways of calculating the mixed derivatives of w, namely 
D qx D qt w(x,t) and D qt D qx w(x,t), the equality of which imposes the constraint 

2V2S 2S2V 2S 2V 

(4 - 26) wi^iiw = D ^ii- D ^w 



A consequence of (|4.25|) is the following companion result to Lemma 14.11 

Lemma 4.3. The Stieltjes function, f, satisfies the q-difference equation 
(4.27) RD qit f = 2Sf + T, 

where T(x, t) is polynomial such that deg x T < deg x S . 

Another compatibility relation is implied by D q t D q x f = D q x D q t f in conjunction with (|4.8[> and (I4.27[) . This 
relation can be stated as 

, Jn „. 2V T 2~S U „ T U 

4.28 — =— = D ax D at — . 

v ' WR RW R q ' W 

When (|4.26p and (I4.28[) are satisfied a companion result to Theorem 14.21 can be stated. 

Theorem 4.4. The matrix Y n is a solution to 

(4.29) D q . t Y n = B n Y n , 



where 

(4.30) B v 



- S -a n $ n 
(R-2t(l-q)S) 9 n - 1 -S-(x-b n -i)^n-l, 



with and ^ n ■polynomials in x specified by 

(4.31a) $„ =R (e n D q>t p n - p n D qtt e n ) + 2Se n p n , 

(4.31b) =a n R n (e n - 1 D q , t p n ~p n -iD q . t e n ) + S(2a n e n -ip n - 1) 
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IS 



PROOF. Our strategy is to adapt the proof of Theorem l4.2l to the D qt operator. Using RD q t f — 2Sf + T and 



p n p 



+ — we have 



RD q t f = RD q t ( rnzi + ±L) =2S 0n - 1+€n +T. 

I'n Pn ' 



This suggests that we define 



(4.32) $„ =R {p n D qtt 4> n -i - (f) n -iD q .tp n ) - 2S(/) n -ip n - Tp n p n , 

=R (e n D q<t p n - p n D q . t e n ) + 2Se n p n . 

Expression (|4.32[) tells us <£>„ is a polynomial in x while (|4.31a[) implies a bound on the degree. 
Using (I4.31a|) and (|4.13[) we arrive at 

a n {pn-l<t>n-\ - Pn<t>n-2)$n = R (PnD q> t4>n-l ~ <t>n-\D q,tPn) ~ 2S(j) n -ip n - Tp n p n . 

By splitting this expression into terms divisible by 4> n -\ and p n , we arrive at an equality that defines given by 

n _ip n * n =0„_i (a„p„_i$„ + (R - 2t(l - q)S)D q ^Pn + Sp n ) , 

=Pn {a n <f> n -2®n + RDq,t<f>n-l ~ S(f> n -l ~ Tp n ) . 

The first line is just a rearrangement of the required ^-difference equation, in t, for p n . The second expression is 
equivalent to 

,t. &n ( Pn~2 ®n RDq,t4>n-l Q Tp n 



<Pn-l <Pn-l 9n-l 

Ra n (j) n -2PnD qt t4>n-l Ra n 4>n-24>n-lD q jPn 2Sa n <t>n-2<t>n-lPn 



Ta n <t>n-2PnPn , RDq^n-l _ g _ Tp r , 



0n-l <t>n-l 9n-l 

=a n R(p n -lD q ,t(t>n-l - 4>n-2D qtt Pn) - S (2a n (f>n-2P n ~ 1) ~ a n Tp n -lPn- 

We remark that this, being a linear combination of polynomials, implies fy n is a polynomial in x. Using (|4.31a[) in 
the first expression for ^ n allows us to write 

Ra n p n -ie n D q , t pn „ „ , 2Sa„p n -ie n p n RD q , t p n 2t(l - q)SD qit p n , a 
W n = Ra n p n -iD q je n H 1 h S, 

Pn Pn Pn Pn 

which upon noting a n p n -\t n = a n p n t n -i — 1 implies (|4.3f bl) . 
The working to date shows 

(R - 2t(l - q)S)D q>t p n = (*„ - S)p n - a„$„p„_ 1 . 

Replacing n by n — 1 in this expression, then using (|2.2j) to express p n -2 in terms of p n and p n -i, establishes that 
p n and p n -i form a column vector solution of (|4.29[) . 

The derivation of the g-difference equation in t for e n /w may also be derived in an analogous manner to the 
proof of Theorem 14.21 so we refrain from the giving the details. □ 
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We note that D 9jt does not necessarily alter the degree in x, hence, the upper bounds for the degrees of $ and 
VP are given by 

(4.33a) deg^, $„ < max(deg x S — 1, deg^, R — 1, 0), 

(4.33b) deg x < max(deg a . S, deg x R, 0). 

Further to this, we may use (12. 6[) to determine coefheients in terms of the aj and &i. 

Equation (l4~29")) may be rewritten in the context of connection preserving deformations to read 

(4.34) Y n (x,qt) = (I - t(l - q)B n {x,t))Y n = B n (x,t)Y n (x,t). 

We use this relation and (|4.18[) to deduce 

P n (x, qt) =(Foo,„(a;, qt)y 1 Y a , n (x, qt), 

=(Yoo, n (x, t))B n (x, ty 1 B n (x, t)Y 0<n (x, t), 

which shows us that the connection is preserved under deformations in t. 

Since Y n satisfies (|2.10|) and (|4.29|) we have a compatibility condition, which follows from a consideration of 

Dq,tYn+l, 

(4.35) M n (x,qt)B n (x,t) + D q>t M n (x,t) = B n+1 (x,t)M n (x,t). 
The first row of (|4.35p is equivalent to 

(4.36a) [R - (1 - q)t(S + * n+1 )] + (1 - q)ta n+1 <5> n+1 = 

X - hn [R-(l- q )t(S + ^ n )}+ {1 - q)tanaA ^ 



"(1 



(4.36b) (1 - q)t{x - 6„)$J + [R-(l- q)t(S + * n+ i)] 



^2- [i? - (1 - g)t(S + *„-i - (x - bn-i)**-!)] 

»n+l 

We have an additional relation 

det Y n a n w a n R 



(4.37) det(J - t(l - q)B n ) = 



det Y n d n w a n (R — 2t(l — q)S) 
A consequence of this relation is the first order recurrence relation in n, given by 

(4.38) a n ((g - l)t (*„ + S) + i?) ((g - l)i (6„_ x - z) + * n _ a + S) + R) 

+ {q - l)H 2 a 2 n $ n - 1 $ n a n - Ra n (2{q - l)St + R) = 0. 

The compatibility condition (|4.19[) is naturally paired with (|4.35[) . Thus, rewriting these read 

(4.39a) A n+1 (x,t)M n (x,t) =M n {qx,t)A n (x,t), 

(4.39b) B n+1 (x,t)M n (x,t) =M n (x,qt)B n (x,t). 
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A further identity which can be grouped with these follows from the compatibility imposed by the requirement that 
D q ,tD q ^ x Y n = D q<x D q<t Y n . One computes 

(4.39c) A n {x, qt)B n (x, t) = B n (qx, i)A n (x, t), 

which is equivalent to (|3.6p . 

5. Deformed little o-Jacobi Polynomials 



The little g-Jacobi polynomials were introduced by Hahn |28j . This family of polynomials possesses the or- 
thogonality relation |44j 

f 1 x a (qxb;q) 00 

Jo (<7^;9)oo 

This ratio of two exponential factors may be scaled and chosen appropriately so that the root and pole is at 0,3 and 
04 respectively. We now adjoin two roots that are proportional to t, ait and a 2 t, to give the deformed weight 

a ( 

(5.1) w(x,t)~ 



( X X . 

ta 2 ' a 4 ' 



In keeping with the notation of [40] . we trust there is no ambiguity between the terms in the three term recursion 
relation, a„, and the roots of the determinant, a;. The deformed polynomials associated with (|5.1[) satisfy 



L{PlPj) — I w(x,t)pi(x,t)p 3 (x,t)dqX = 5ij. 

Js 

The set S, also called the support of the weight, may begin and end at distinct roots of w(x,t). These include 
03, a\t and 0. Choosing 0,3 and ait and using (|4.5[) allows the moments to be expressed in terms of Heine's basic 
hypergeometric function, 

. ~a+k fx X . _ \ 

f«» t x {^^>9) 

Mfc = / 7 x d q X, 

V ta 2 ' a 4 ' 



[q aity+*+\i- q )(^q^, q - q ) 



04 n a+k+l 

a 3 > y 



( n o+k+l qai. 
\ q ' a 2 ' 



q ait 

9; 

04 



\ * / OO J. I CL\ 1 1 



, 903 

' a 2 t 



/ 00 v 

This allows us to use (|2.5|) to express a n and b n in terms of determinants of basic hypergeometric functions. 
This weight (|5.1|) satisfies the equation 

_ / a 2 a 4 (a; - ait)(a; - a 3 ) - q a aia 3 (x - a 2 t)(x - a 4 ) \ 
Dq,x w \ x i<>) — I . . , . . . I W(X,t). 

\ a 2 a4(x — ait)(x - a 3 )x(l - q) J 

A comparison with (14.71) reveals that the spectral data polynomials are 

W =a 2 a^{x — ait)(x — a^)x{l — q), 

2V =a 2 o,i{x — ait)(x — 03) — q° aia 3 (x — a 2 t)(x — 04). 
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Recalling Theorem 14. 2 [ it follows that the poles of the linear q- difference equation in x satisfied by these polynomials 
is determined by the polynomial 

(5.2) W - 2x(l - q)V = q a a 1 a 3 (x - a 2 t)(x - a 4 )x(l - q). 
In the t direction, w satisfies the equation 

( ai{x - qa 2 t) - a 2 {x - qait)\ 

Dq, t w(x,t) = — r W{x,t). 

\ ax{x - qa 2 t)t{l - q) J 

Comparing this expression with (|4.25p shows 

R(x,t) =a\{x — qa 2 t)t{\ — q), 
2S(x, t) =a\{x — qa 2 t) — a 2 (x — qa.it). 

The appropriate poles of the linear q-difference equation in t satisfied by these polynomials is therefore determined 
by the polynomial 

(5.3) R - 2t{l - q)S = i(l - q)(x - q ai t). 

We remark these explicit forms for W, V, R and S satisfy (|4.26[) as they must. 

5.1. Linear problem. Since we have an upper bound for deg^. 0„, deg^ r2 n , deg^, $„ and deg^, * n from (|4.f 61) 
and (|4.33j) . we parameterize 6„, O n , $„ and ^ n by 

(5.4a) e„ = O , « + i,„x, 

(5.4b) 0„ = UJo.n + Ul,n& + ^2,nX 2 , 

(5.4c) $„ = cf)o, n , 

(5.4d) = ipQ. n + ipi in x. 

This bounds the degree of the relevant polynomial component of the linear q-difference equations in x and t. Hence 
the linear q-diffcrcncc equations satisfied by the polynomials may be written in the form (|4.17[) and (|4.34l) where 

A ,„ + A, n x + A 2 „x 2 



(5.5) A n = I — x(l — q)A n 

(5.6) B n =I-t(l- q)B n 



(x ~ a 2 t)(x — 04) 

Bo,n + Bl n X 



(x — qa\t) ' 

for some set of Ai j7l and Bi^ n . According to (|4.22l) and (14.37[) , the determinants of these matrices are 

a 2 a^{x — a\t){x — an) 



det A, 
det B n = 



aia3q cr (x — a 2 t)(x — 0,4) ' 
aia n (x - a 2 qt) 



a 2 (x - aiqt) a n 

At this point, the associated linear q-difiercnce equation satisfied by the orthogonal polynomials is one in which 
the coefficient matrix, A n , is rational rather than polynomial. To relate this formulation to the classical theory 
of Birkhoff [10] . or more precisely, Jimbo and Sakai [40], we require a gauge transformation that will relate the 
linear q-difference equation in which the coefficient matrix is rational to a linear q-difference equation in which the 
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coefficient matrix is polynomial. By considering the associated q-difference equation for Y* = Z n Y n , we note that 
Y* satisfies the trio of equations 



(5.7a) Y^qx, t) = (Z n (I - x(l - q)A n )Z^) Y 7 l = A^Y*, 

(5.7b) Y:(x, qt) = (Z n {l - t(l - q)B n )Z- 1 ) Y* = B*Y*, 

(5.7c) Y; i+1 (x, t) - (Z^MnZ- 1 ) Y* = M*Y*. 

By letting Z n to be proportional to appropriate g-exponential factors allows A* to be polynomial. We may also 
choose Z n carefully so that A* possesses some desirable properties, such as certain asymptotic characteristics in x 
and/or t, and doing so makes the correspondence to the work of Jimbo and Sakai [40] more apparent. Specifically, 
by choosing 

/ a 2 a 4 q-" \ 



Z n {x, t) 



( x x 
02 1 ' 0.6 

we have that Y* satisfies the g-difference equations 





e q,qa,2 WTn 

7 "- 1 



\ e q,qax (t) } 



(5.8a) Y^qx, t) =(A* „ + A* l n x + A^x^ = A;r„*, 

(5.8b) Y:(x, qt) = g ( fl 5,n + giV) Y . = b , y: 

(x — qa\t)\x — qa 2 t) 

The corresponding determinants are given by (14.22[) . (|4.37l) and (14 .6|) 

(5.9a) det(A*) =a\a 2 aza^q a (x — a\t){x — a,2t)(x — a^)(x — 04), 

(5.9b) det(B*) =- -. 

y ' y n> (x - qta l )(x - qta 2 ) 

It will transpire that the form of the coefficient matrices of (|5.8aj) and (|5.8bj) is well suited for the purpose of 
parameterizing the linear problem satisfied by the orthogonal polynomials. 

Although (|5.8al) specifies a 2 x 2 linear (/-difference system in which the determinant of coefficient matrix, given 
by (|5.9a|) . has roots that coincide with those found in |40) . we require two additional properties; firstly that A* 2 n is 
a constant diagonal matrix and secondly, that Aq n is semisimple with eigenvalues proportional to t. 
An asymptotic expansion of Sl„ and 0„ around x = 00 reveals 

1204 — q a (2q n — 1)0103 



U 2 n 



— -q n+ ° ai a 3l 



giving 

(5.10) Al n = 

where 




Ki = q n+a aia 3 , k 2 = a 2 a 4 q 



This shows that the linear problem possesses the first required property. 
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To show that n has eigenvalues that are proportional to t, we first write Aq n as 

(KiK2t + a^a^a^t qK2W n 9 . n 
o u ^,n 
~ O n -l"0,n-l ~ w ,n-l 
qK2W n 2 

where we have used the notation 

(5.12) Wn = e ^iW(^-fi) . 

qe q ,qa 2 {t)X 

We let Ai,„ and A 2 ,„ be the eigenvalues of Aq n . Utilizing (|5.9[) and that det(j4* (0, t j) — \i, n \2, n gives us 

Al,nA2, n = KlK2aid2«3C4i 2 ; 

revealing that Ai, n A2, n = q 2 Xi. n X 2 , n . Adding Bq n to both sides of the residue of (|4.39cp at x = 0, namely the 
relation 

A),n^o,n = qB Qn A 0n , 

and then taking determinants shows 

det(/ + iS,J =det(J + gAS ! „), 

revealing 

1 + Ai, n + A2,n + Ai, n A2, n = 1 + ?Ai jn + gA2, n + (? 2 Ai in A2,n- 

This shows that Ai i „A 2in oc t 2 and Ai, n + A2. n oc t, hence Ai jn and A2. n are proportional to t. 

A further property of Ai in and A 2 .„, that will be useful later on, is their independence of n. The independence 
of kik 2 's on n indicates that Ai )n A2, n is independent of n. However the trace of A^ n is 

K\K 2 t ( 1 H -J — (u!0,n + WOn-l + b n -l6o n — l) = Al n + ^2,n, 

\ <f ) 

which indicates that the constant coefficient of (|4.20[) may be expressed in terms of Ai.„ and A2, n as 

Ai, n + A 2j „ — Ai )n +i — A 2j „+i = 0. 
This proves Ai,„ + A 2 ,„ is independent of n, hence Ai in and A2, n are independent of n. We may now write 

(5.13) {Ai,„,A 2 ,„} = {0^,02*}, 

where 9\ and 02 are constant in t and n. These eigenvalues are not free, with an implicit dependence on the a^s 
and Kj's and the support chosen. 

The additional properties mean that (|5.8ap can be cast in a form equivalent to the linear problem studied in 
[40] . Technical achievements in [40] are to identify the parameterization of the linear problem which leads to the 
q-Pvi system. The present orthogonal polynomial setting allows us to perform these steps in a more detailed, and 
perhaps more systematic manner. 
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5.2. Orthogonal polynomial parameterization. Our pathway toward the parameterization of the problem 
is to make use of the orthogonal polynomial variables. Parameterizations of this sort can be found in previous works 
® [HI 1231 124j . However, these works do not provide a systematic way to link up with co-ordinates that specify 
Painleve systems. 

To begin, using the expansions (|4.11l) of Q n and 0„, we find 



(5.14a) 
(5.14b) 



(5.14c) 



(1 — q)ni " 

wi,n = hi + ki (a 2 t + di) - '" 1 "' z r) v ~^ ' — - ^a 2 a 4 (ait + a 3 ) . 



i=0 



K\K% (a 2 t + a 4 ) 1 



'0, n — - 



Ki (a 2 qt + a 4 q-q Yh=q h + Yn=o & «) K 2 (aiqt + a 3 q + q ^=o b i - E"=o b > 



/n-l 



- m a? 



a 4 Ki I ^ h I - — — 



3 

I'n-l 



n— 1 n— 2 n— 1 

Kl \ -s 2 K l \ \ - , , K\ 



q 

n—1 n—1 



i=0 



\i— 1 / z—1 i— i=i+l i— 

2 



— ( ^ 6,: ) H - a 2 a 4 tK 1 



tn-1 



K 2 2 



a\a 2 a 3 a 4 t tK\K 2 



The expansions (14.311) of $ n and give 
(5.15a) 



lpi,n =77 («1 + 02) 

2 7n 



(5.15b) 
(5.15c) 



7" 



a 1 a 2 qt 1 



fo.n 



017^ - a 2 7^ 



7rx7n 

This specifies a parameterization of the linear problem for 1^ in terms of orthogonal polynomial variables. We use 
the notation 

n-l 

(5.16) r„ = ^& i5 

i=0 

which is proportional to the coefficient of x n ~ x in p n . By combining (|5.14l) . (|5.4I) . (14.101) and (|5.8al) . 



/ (q - l)«ir n _ 



(5.17) 



A* — 

A l,n — 



K\ (a 2 t + a 4 ) 



K 2 W„ 



\ 



a n (qui - k 2 ) (qn 2 - m) (q - l)K 2 r r , 



V 



q 2 K 2 w n 



k 2 (ait + a 3 ) 



We make use of the relations 



traced n =0 4 t + 8 2 t, 



det A 



Q.n —^1^2* , 



which allows (|5.1ip to be simplified to 

/ t (aia 2 a 3 a 4 + KiK 2 ) 



(5.18) 



x 0,n 



V 



(qKx - k 2 ) al9 Q . n -i 



qn 2 w n 



qK 2 W n 0o,n 

qui — k 2 
t (a\a 2 a 3 a 4 - 2 (6% + 8 2 ) + kik 2 ) 
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where 

{a\a2a^a4t — 2t6\ + tniK2 — 2wo,n) (01020304^ — 2t02 + tn\K2 — 2wo,«) 
We simplify the expression for Aq n by introducing the variable r„ so that we may write Aq n as 



qui - k 2 



Oit + r n 

{:>A<}) A ln=\ ( gKl - K2 )r tl (r n + M-g2t) 

qK2W n 9 , n 

In relating the coefficient of x 2 in det j4* with (|5.9a[) , we express r n as 

_ {q - 1)kik 2 (ait - a 2 t + a 3 - 04) T n _ (1 - q) 2 niK 2 T 2 
q(Ki-K 2 ) q 2 {Ki-K 2 ) 

a 2 n (qui - n 2 ) (qn 2 - K\) t (6 2 Ki + 9 1 k 2 - aia 3 K 2 Ki - a 2 a i K 2 Ki) 
q 2 (k\ - k 2 ) hi - k 2 

Equating (|5.19j) with (|5.I8|) gives an alternate representation of aJo,n to that of (|5.14c[) . The equations (|5.19j) . (|5 . 1 7|) 
and (|5.10|) are explicit parameterizations of the linear problem using orthogonal polynomial variables combined 
with knowledge of the structures (|5.8a|) and (|4.I0p . 

We now turn our attention to the parameterization of the linear problem, (|5.8b[) . involving B* n . First, upon 
recalling (I4.30[) . it follows from the large x expansion of and 1 $ n , as implied by (I5.15|) . that 

B* hn = -tl. 

Direct substitution of the values of ^> n and <&„ from (|5.I5[) gives 



(5.20) B* = t 



T„ - r„ + a x qt 
(qui - k 2 ) (w n a 2 n - a^Wr, 



qui - k 2 
r n - f „ + qa 2 t 



\ qK 2 w n w n J 

This gives us enough information to deduce the evolution of the variables j 2 , a\ and T„, which completes the 
parameterization of the linear problem in terms of the orthogonal polynomial variables. 

Use will be made of (|5.fOI) . (|5.17p . (|5.19[) and (|5.20[) as we now proceed to make the correspondence between 
the above discrete dynamical system and <?-Pvi by making a correspondence between the linear systems. 

5.3. Jimbo-Sakai parameterization. Our primary task is to find expressions for w n , y n and z n in terms of 
7^ , a 2 n and b n and vise versa. We have chosen w n in the previous parameterization, as it is related to j 2 via (|5.I2|) , 
to be the variable that reflects the gauge freedom in both parameterizations of the linear problem. In keeping with 
earlier remarks, we deduce 

(5.21) e , n = y^~* 2 \ 

q 

and define variables z\ and z 2 according to 
(5-22) A* n {y n ,t) = 

Evaluating the determinant at x — y n reveals 

(5.23) ziz 2 = (y n - ait)(Vn ~ a 2 t)(y n - a 3 )(y n - a 4 ). 
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We factorize this into the factors 

(5.24a) Zl Jlfa-taiXlfa-taa) 

qniz n 

(5.24b) z 2 =(y n - a 3 ){y n - ai)qniz n . 

The benefit of this particular factorization reveals itself in the proof of Theorem 15. II It follows from (|5.10j) . (|5.11|) 
and (|5.19[) that z\ and z 2 may be expressed in terms of a\ and T n via the expressions 

ft . ok \ KiQo,n ((? ~ l)r n - q (a 2 t + a 4 )) ? 2 «i0o,n . . , a 

(5.25a) KiZi = ; V -. — +r n + t9i, 

n 2 -qni \K2-qKi)- i 

/, „, u K 2 o , n ({q- l)T n + aiqt + a s q) q 2 K 2 6 2 

(5.25b) k 2 z 2 = h - 2 -r n + t9 2 , 

n 2 — qn\ ( K2 — (jf/tjj 

which specifies z n . To be consistent with (15.221) . the matrix in (I5.8aj) permits the parameterization [40] 



(5.26) A* r: 



' Ki{{x - y n )(x - a) + zi) K2W n (x~y n ) 



Kifax + d) ,, . . . 

K2{{x - y n ){x - p) + z 2 ) 



w„ 

where a, (3, 7 and 6 are to be determined. Comparing the upper left entry of (|5.26l) with (I5.10p . (15.171) and (|5 . 19[) 

shows 

q (a 2 t + a 4 - y n - a) 



(5.27) r„ 



These substituted into (|5.25b[) reveal 



q-1 

=Kiy n a + kizi - t6i. 



(5.28) a = [ — {{9i + 2 )t - k\Zi — k 2 z 2 ) - n 2 ((ai + a 2 )t + a 3 + a 4 - 2y n ) 

K1 — K2 \ Vn 

Conversely, comparing coefficients of the lower- left entry of (|5.26l) with (|5.10[) . (|5 . 1 T[) and (|5 . 19[) gives 

q(a 1 t + a 3 - y n - 0) 



(5.29) r 



r„=- K 2 y n /3 - k 2 z 2 + t9 2 . 



These substituted into (|5.25al) show 



(5.30) (3 = [ ((9i + 6 2 )t - k\z\ - n 2 z 2 ) + Kx((ai + a 2 )t + 03 + a 4 — 2y n ) 

«1 - «2 V Vn 

The strategy to be used to specify 7 and 5 makes use of (|5.9a|) . By equating the coefficient of x 2 of det A* n from 
([5T2^|) with ([Oa]) . we have 

(5.31) 7 = zi + z 2 + (y n + a)(y n + j3) + (a + j3)y n - aia 2 t 2 - {ax + a 2 )(a 3 + a 4 )t - a 3 a 4 . 
A comparison between the coefficient of x in (|5.26[) with that of (|5.9a[) shows 

(5.32) 5 = — (aia 2 a 3 a 4 t 2 - (ay n + Zi)(/3y n + z 2 )) . 

Vn 

This concludes our task of parameterizing the linear problem associated with the orthogonal polynomial system 
with the weight (|5.ip and its correspondence with the parameterization of |40j . After completing the task of 
parameterization of a, (3, 7 and S, Jimbo and Sakai proceeded to give the coupled equations referred to as q-Pyi- 
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However few details were given there. We shall provide details by making use of a structured form of the B matrix 
following from the orthogonal polynomial viewpoint. The structured form of the B matrix follows by using the 
substitutions of ([537]) and (fQ9|) in (|O0|) . giving 



(5.33) 



D* 

B 0,n - 



qt 2 (ai + a 2 - D q j(y n + a)) 



qtn 2 (w n - w n ) \ 



qui — k 2 



\ 



qt 2 (oi + a 2 - D qt t(y n + /?)) 



(«i - qn 2 )w 

In addition to (|5 . 33[) a crucial ingredient in our derivation of q-Pyi are the compatibility conditions (j4.39[) . After 
making the transformation Y* — Z n Y n these latter conditions read 



(5.34a) 
(5.34b) 
(5.34c) 



A^(x,qt)B*(x,t) = B*(qx,t)An(x,t), 
A* n+1 (x,t)M*(x,t) = M*(qx,t)A* n (x,t), 
B* +1 (x, t)M* n (as, t) =M*(x, qt)B* n (x, t) . 



By evaluating the residue of (|5.34ap at x = ait, a 2 t, qciit, qa 2 t, we obtain the expressions 

(5.35a) (qaitBt >n + B^A^t, t) =0, 

(5.35b) {qa 2 tB* l n + Bl n )A* n {a 2 t, t) =0, 



(5.36a) A* n (q ai t, qt){qa x tBl n + B* n ) =0, 

(5.36b) A* n (qa 2 t, qt)(qa 2 tB* n + B*J =0. 

By looking at the residue of (|5.34a[) at x = 0, we obtain the additional relation 

(5-37) A 0n B Qn = qB Qn A Qn . 

Theorem 5.1 ([40]). The compatibility condition, (|5.34al) . is equivalent to the evolution equations for y n , z n and 

w n specified by 

{qnxZn - 1) 



(5.38a) w n 
(5.38b) 

(5.38c) y n y n 



K2Z n - 1 

(y n ~ a\t){y n - ta 2 ) 
qKiK 2 {y n - a 3 )(y„ - o 4 ) ' 
q(0\z n - taia 2 ){8 2 z n - ta\a 2 ) 



aia 2 {qniz n - l){n 2 z n - 1) 

Proof. For brevity, we let the parameterization of B^ of (|5.33p be given by Bq = (&jj)jj=i,2- The upper 
right entries of compatibility condition (|5.36a[) and (|5.36b|) read 

K 2 w n (y n - qtai)(b 22 - qait 2 ) = Ki&i 2 ((gtai - y n ){qtai - a) + Zi), 
K 2 w n {y n - qta 2 )(b 22 - qa 2 t 2 ) = K 1 bi 2 {{qta 2 - y„){qta 2 - a) + zi), 

which gives us an expression for 612 and 622- We only require 612, which is given by 

tn 2 w n (y n - qait)(y n - qa 2 t) 

012 — 



Ki({y n - qa\t){y n - qa 2 t) - zi) 
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Equating this with the upper right element of (|5.33l) gives 

q(w n - w n ) w n (y n - qait){y n - qa 2 t) 



qui ~ n 2 ki(zi - (y n - qait)(y n - qa 2 t)) ' 

This evolution equation is simplified using the particular factorization (|5.23[) . The structure of the right hand side 
of the above relation justifies, a posteri, the factorization (|5 . 23|) . The particular form of (|5.24|) means the evolution 
of w n is equivalent to (|5.38a|) . 

The upper right entries of compatibility condition (|5.35a[) and (|5.35bj) read 

n 2 w n (y n - M)On - qa\i 2 ) = K 2 b 12 {{ta 1 - j/ n )(iffli - 0) + z 2 ), 
K 2 w n (y n - a 2 t)(bn - qa 2 t 2 ) = K 2 bi 2 {{ta 2 - y n )(ta 2 - /3) + z 2 ), 

which we solve in terms of bn and 6 12 to give 

qtw n (y n - a\t)(y n - a 2 t) 



(5.39) 
(5.40) 
We deduce 



b\ 2 = 
b n = 



{Vn - ait)(y n - a 2 t) - z 2 ' 

qt {z 2 {y n - (ai + a 2 ) t) + j3 {ait - y n ) (a 2 t - y n )) 
(ait - y n ) (a 2 t - y n ) - z 2 

K 2 (w n ~w n ) w n (y n - aii)(y„ - a 2 t) 



qni-K 2 [y n ~ ait)(y n - a 2 t) - z 2 
which is equivalent to (|5.38b[) knowing (|5.38al) . Comparing (|5.40[> with (|5.33l) yields 

■ / . , rs , xx z 2 (y n - (ai + a 2 ) t) + fi (ait - y n ) (a 2 t - y n ) 

t [a\ +a 2 + D q , t (y n + a)) = — — — , 

(ait - y n ) (a 2 t - y n ) - z 2 

which is equivalent to (|5.38c[) knowing (|5.38a|) and (|5.38bl) . or the particular Riccati solutions 

. _ gy n (l - K 2 z n ) 

Vn — 1 x , 

1 - qn 2 z n 

the latter not being satisfied in general. The derivation of the evolution equations is complete. 



Full correspondence with the Jimbo and Sakai form is obtained by letting 



(5.41) 

where (|5.38[) become 



«5 



aia 2 



ae 



aia 2 



a-j = 



qui 



1 

K 2 ' 



y n y n 



under conditions that 



as given in [40] . 



a 7 a 8 (y n - a\t)(y n - ta 2 ) 
(y n - a 3 )(y n - a 4 ) 

a3a±(z n - a<jt)('z n - a 6 t) 
(z n — a 7 )(z n - a s ) 



a5«6 _ qaia 2 
a 7 a 8 a^a4 



□ 
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We now return to the orthogonal polynomial context for these results. In addition to (|5.7ap the three term 
recursion relation, (|2.10[) . in the orthogonal polynomial context gives us another linear problem. The representation 
of M* following from flUTO} and (jlTfc]) is 

( x - b n n 2 w n 
A qKl Q K2 
qn 2 w n 

This can be used to express the orthogonal polynomial quantity b n in terms of the natural variables. Considering 
the coefficient of x 2 and x in the upper left and right entries of (|5.34bj) respectively results in the expression 

(5.43) bn = g(g : ia " K2/3) . 

q-Ki - n 2 

For the orthogonal polynomial quantity a comparison of the lower left component of A\ n given by (|5 . 1 T[) and 
(fS36|) shows 

(5.44) al = *«™ 

{qui - K 2 )(qn 2 - Ki) 

One important consequence from this perspective is that the natural variables may be expressed in terms of 
determinants of the moments. Using (|5. 12|) and (|2.5a[) we have 

(5.45) w n = e,Wflte-g"i)*n-K 
Using (|5.14bl) and (I5T2T]) gives 



, . qni(a 2 t + 04) - K 2 (ait + 03) Ki — K 2 S n l^i ~f S„ +1 

(5-46) y„ = + -r --r . 

qui - k 2 qui - n 2 A„ qui - n 2 A n+ i 

The simplest dctcrminantal form for z n comes from the substitution of (|5.45p into the inversion of (|5.38a[) . which 
reveals 

aiA„ + iA n — a2A„A„ + i 

(5.47) 2 



A 

n+ 1 



These may correspond to known determinantal solutions, such as the Casorati determinants of Sakai [56], although 
we are yet to investigate this point. 

5.4. Backlund transformations. The linear problem equivalent to the orthogonal polynomials three term 
recursion, (|5.7c[) , may be expressed in terms of the natural variables appearing in (|5.43[) . Substitution of (|5.43[) into 
HElg) gives 

/ q 2 Ki(x — a) + K 2 (qj3 — x) n 2 w n 

M* = 



V 



q 3 Ki - qn 2 qn\ - k 2 

k 2 - qui 



qn 2 w n 

In the context of orthogonal polynomial theory the system of equations describing the evolution of this system in the 
n direction are known the Laguerre-Freud equations. Moreover, these very recurrence relations in the transformation 
n — > n— 1 and n — > n+ 1 represent elements in the group of Backlund transformations. Since the group of Backlund 
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(5.48) 



Vn-l Zn-1 



transformations are of affine Weyl type, the Laguerre-Freud equations are equivalent to a translational component 
of the extended affine Weyl group of type Z)g . We represent the n n — 1 translation as 

ai a 2 a 3 a 4 a\ a 2 a 3 a 4 

a 5 a 6 a 7 a 8 J a 5 a 6 qa 7 2a 

The derivation of its explicit form relies on (|5.34b|) . The lower right entry of (|5.34b| shifted n — > n — 1 at x — y n -i 
yields the relation 

(5.49) y n _ x = — . 

7 

By evaluating the upper right entry of (|5.34b[) shifted byn->n-lati = Vn-i we obtain 

(Vn-l - Vn)(Vn-l - a) + Zi 



(5.50) 



Zn-l = 



qn 2 (flA - Vn-l) (Vn-l - a 3) 



Finally, using (|5.12l) to find w n -\/w n reveals 



w n -i 



al{qni - k 2 ) ' 

which expresses y n -\, Zn—i an d Wn-i in terms of y n , z n and w n . 

A more canonical transformation from the orthogonal polynomial perspective is the transformation correspond- 
ing to the shift n — > n + 1, which is represented by 

ai a 2 a 3 a 4 I a i a 2 a 3 «4 

: y n z„ V -> < : y„+i z n+x 

a 5 a 6 a 7 a 8 J a 5 a 6 2Z <ja 8 

Another viewpoint is that this shift is a ^-difference analogue of a Schlesinger transformation of the linear system, 
which induces a Backhand transformation of the Painleve eauation |38j . The Schlesinger transformation is induced 
by multiplication on the left by a rational matrix, this rational matrix coincides with M n (x) for this particular 
solution of the linear system. 



Theorem 5.2. The shift {y n ,z n ) — > (y n +i, z n+ i) is given by 



(5.51) 
(5.52) 



z n +i 



K 2 Z n [y n {a\t - y n ) + C«] [y n {a 2 t - y n ) + Cn] 

g 2 Kl |/C 2 y„Z„ (03 ~ 2/n) + Cn] [«2yn2n (04 - J/„) + Cn 



_ ^2j/n (1 ~ K2-g K ) 

q A n\ (1 - q^KiZ n+ i) 



Cn - {y n ~ ati)(y n - a 2 t) + 



z n {q0\t - n 2 a\a 2 t 2 ) 
1 - n 2 z„ 



K 2 y n z n (a 3 -y n ) + („ 



Cn - (y n - ait)(y n - a 2 t) + 



z n (qd 2 t - K 2 aia 2 t 2 ) 

1 - K 2 Z n 



^2ynZn (04 - y n ) + Cn 



where 



(k 2 - q «i)Cn = «2(2/ n - axi){y n - a 2 t) - q KiK 2 {y n - a 3 )(y n - a±)z, 



K 2 z n (tQi - qK 1 a 3 a 4 )(t6 2 - q^a^) 
1 — K 2 z n ma 3 a4 
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PROOF. Using (|5.34b[) . we note that an alternate way of writing A n+1 is given by 

(5.53) A n+1 (x,t)=M*(qx,t)A n (x,t)(M*(x,t)y 1 . 

Using (|5.26p to represent the top row, and the right hand side of (|5.53[) to express the bottom row, we have 

(qni((x - y n +i)(x - a) + zx) q^ 1 K 2 w n+1 (x - y n +i) 
(k 2 - qni) 2 (x - y n ) , , _ x v . — x \ _] 
(x - y n ) [On [Kx - q L n 2 )-ma + q xk 2 ) + z X Ki 
qK 2 W n 

where Zi and a denotes Zi and a at n + 1. The determinant of A n+ i at x = ait is zero. However, using this 
representation of A* l+1 , the top row is divisible by {y n +i — o-it) and the bottom row is divisible by (y n — ait). This 
also applies to the case for x = a 2 t, hence by equating the determinant of A n+i with zero gives two expressions for 

W n +1 

{K KA\ f yn+X-Oit 2 \ ( yn- Ojt . , \ (k 2 - g/tl ) 2 W n+ i 

(5.54) q Ki (ait - a) ) [ (qui - k 2 ) b n + quia - aitn 2 = , 

V z n+l " ) \ Z n J W n 

fKKK s (y n+ i-ait 2 -AfVn-axt , s, . . \ («2 ~ q^x) 2 w n+ i 

(5.55) q Ki (a 2 t - a) (qui - k 2 ) b n + quia - a 2 tn 2 = . 

V. Zn+X J \ Z n J W n 

In a similar manner, we consider the matrix representation of A* n+1 given by 

I (x - y n ) (ni(qx - b n ) + q^ 1 K 2 (b n - qP)) + K2Z 2 q^ 1 n 2 w n+ i(x - y n +x) 

A n+i = (n 2 -qni) 2 (x-y n ) n 2 - 

— {{x - y n +x){x - P) + z 2 ) 

\ qn 2 w n q 

which has been obtained by using the left hand side of ()5.53|) to represent the left column of A n+i and the right 

hand side of (|5.53j) to represent the right column of A n+ i. The left and right columns are divisible by y n — a% and 

y n +i — CI3 respectively at x — a%. This applies also in the case of x = 0,4. Hence equating the determinant of this 

representation of A n+i at x — 03 and x = 04 with zero gives two additional equations for w„+i 

(5.56) - z„(a 4 — y n ) H 1 z n +x (ax - y n +i) H 5 = ~? 

V q \qni K 2 J qK\ n 2 J \ q z n x / q z \qm n 2 



1 


K2) 


_ JL. 


\- — 


qKi 




qni 


K 2 


1 




_ A_ 


aA 

V — 


qui 




qKi 


K 2 



, r , , , . b n ( 1 \\ /3 a 4 \ / a 4 -/A 1 ( 1 1 \ iu„+i 

(5.57) - z„(a 3 -y„)H 1 z n +x (a 3 - y™+i) H ^ = ~2 \ ■ 

q \qni K 2 J qui k 2 J \ q z K X J q z \qni k 2 J w n 

Equating the coefficient of x in the upper right entry of (|5.34bj) with zero reveals 

. , (qKi - K 2 ) b n + q (qniy n + k 2 P) 

a + y n +i = 2 1 

q z Ki 

which reduces (|5. 54115. 57)) to expressions for w n+ i that are all of degree one in y n +i and The compatibility 

between (|5.54[) and (|5.55|) is equivalent to 

. 2 \<2 -a a 3 a 4 KiK 2 (n 2 z n - q 2 Kiz n+ i) (aia 2 t - q9iz n ) (aia 2 t - q9 2 z n ) 

(5.58) y n +iz„ [k 2 -q KiJ (q Kiz n +i - 1) = t-t 7 -7 

^1 6*22/71 (K2Z n - 1) 

- g 2 Ki (z n - Zn+l) ((03 + 04) K 2 Z n - (01 + a 2 ) t) 

+ g 2 y n Kl («22n - 1) (Zn - ^n+l) ■ 

as is that of (15.561) and (|5.57[) . Substituting f|5 -58[) into the equation resulting from the comparison of f|5 . 54[) and 
(|5.56[) the yields (|5.51l) . To obtain (|5.52[> . we substitute the expressions for z n+ i, given by (|5.5ip . into the right 
hand side of (1535)) . □ 
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As a preliminary check of the recurrence relations, we may consider the special case in which the support is 
chosen to be between and qa±t. In this special case the moments are 

{qai tr^(l- q )(^^, q ; q ) 
Mfc - 7 \ 

We explicitly compute the eigenvalues of Aq to be 

9i = <f '01020304, 9 2 — 01020304. 
Substituting these values of /ik into (|5.46|) . (|5.47|) . (|2.3|) and (|2.4[) for the n = case gives us the seed solution 



Vo = 
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As an illustration of the computation content of the recurrence relations and as a check on their veracity, we 
may compare numerical values of y n +i and z n+ \ using (|5.46|) , (|5.47l) , (|2.3[) and (I2.4[) found by using (|5.52l) and 
(|5.51|) from y n and z„ for generic values of the parameters, t and small values of n. Numerical evidence has been 
obtained to verify that (yi, zi), found using (|5.46[) . (|5.47[) . (|2.3j) and (|2.4[) . coincides with the values of (yi, Zi) found 
by using (|5.52|) and (|5.51|) from the values of (yo,zo) and (|5.46[> . (I5.47p , (|2.3p and (|2.4|) . In a similar manner, we 
were also able to test the relationship between (yi,zi) and (y 2 ,z 2 ) using (|5.52| and (|5.51|) compared with values 
obtained by using (f5T46|) . (|Q7|| . ([2~3]) and (|2~4"1) . 

We remark that the evolution n — > n + 1 of the linear system corresponding to a deformed version of the Pastro 
weight supported on the unit circle, which is the circular analogue of the little g-Jacobi weight, has recently been 
obtained by Biane 9 . The structure of the iterations in n should have a similar structure to other translational 
components of the affine Weyl group, such as the translational component that coincides with the evolution of 
g-Pvi- This multiplicative structure is of (|5 . 52[) and (|5.51[) is similar to the Backhand transformation of Biane [9]. 
In the work of Biane [9], the Backlund transformation, representing the shift n — > n + 1, simultaneously changes 
one of the eigenvalues of Aq n and A 2 n , whereas in our transformation, the eigenvalues of Aq n are independent of 
n. The little g- Jacobi case has also been studied in [27] , although in a truncated way. It is clear from this work that 
the authors have treated a specialized, in the sense that t is fixed by the parameters, and a degenerate case, whereby 
the parameters are related by 0104 = 0203, and consequently have recovered elementary function expressions for 
the three-term recurrence coefficients. 
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